Quantum interference in the classically forbidden region: A parametric oscillator 
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We study tunneling between period-2 states of a parametrically modulated oscillator. The tun- 
neling matrix element is shown to oscillate with the varying frequency of the modulating field. The 
effect is due to spatial oscillations of the wave function and the related interference in the classically 
forbidden region. The oscillations emerge already in the ground state of the oscillator Hamiltonian 
in the rotating frame. 
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Nonlinear micro- and mesoscopic vibrational systems 
have attracted much interest in recent years. In such sys- 
tems damping is often weak, and even a comparatively 
small resonant field can lead to bistability, i.e., to coex- 
istence of forced vibrations with different phases and/or 
amplitudes. Quantum and classical fluctuations cause 
transitions between coexisting vibrational states. The 
transitions are not described by the conventional theory 
of metastable decay, because the states are periodic in 
time and the systems lack detailed balance. Experimen- 
tally, classical transition rates have been studied for such 
diverse vibrational systems as modulated trapped elec- 
trons Josephsoni unctions Q, nano- and microme- 
chanical oscillators 0, Q, Q , and trapped atoms Q , and 
the results are in agreement with theory 0, Q . 

Currently much experimental effort is being put into 
reaching the quantum regime 0, [Io| . In this regime 
tunneling between coexisting classically stable periodic 
states should become important, for weak dissipation. It 
was first studied for a resonantly driven oscillator, where 
a semiclassical analysis [llj made it possible to find the 
tunneling exponent in a broad parameter range 

Tunneling is particularly interesting for a parametri- 
cally modulated oscillator. Here, the coexisting classical 
periodic states have period 2rp , where Tp is the modu- 
lation period. Such period-2 states are identical except 
that the vibrations are shifted in phase by n. There- 
fore the corresponding quantum states (Floquct states) 
are degenerate. Tunneling should lift this degeneracy, 
as for a particle in a symmetric static double-well po- 
tential. Earlier the tunneling matrix element was found 
fl3| for modulation at exactly twice the oscillator eigen- 
frequency loq. Recently the tunneling exponent was ob- 
tained in a general case where the modulation frequency 
ujp = 2-k/tf is close to 2lu [3] 

In this paper we show that tunneling between period-2 
states of a parametrically modulated oscillator displays 
unexpected features. We find that the tunneling matrix 
element oscillates with varying u)p — 2luq, periodically 
passing through zero. These oscillations are accompanied 
by, and are due to spatial oscillations of the wave function 



in the classically forbidden region. 




FIG. 1: (Color online) The scaled effective Hamiltonian of the 
oscillator in the rotating frame g(Q,P), Eq. ([SJl, for /j, = 0.5. 
The minima of g(Q, P) correspond to the period-2 vibrations. 
The eigenvalues of g give scaled oscillator quasienergies. 



For resonant modulation, \u)p — 2w | "C wy, and for a 
small amplitude of the modulating field F the oscillator 
dynamics is well described by the rotating wave approx- 
imation (RWA) [IH. The scaled RWA Hamiltonian g 
as a function of the oscillator coordinate Q and momen- 
tum P in the rotating frame is independent of time. In 
a broad parameter range it has a symmetric double-well 
form shown in Fig.[TJ The minima correspond to the clas- 
sical period-2 states, in the presence of weak dissipation. 
Respectively, of utmost interest are tunneling transitions 
between the lowest single-well quantum states of g. 

A simple model of a nonlinear oscillator that describes 
many experimental systems, cf. Refs. [l|, [H, 0, HI HJ @, is a 
Duffing oscillator. The Hamiltonian of a parametrically 
modulated Duffing oscillator has the form 

H = \p 2 + \ H + Fcosuopt) q 2 + i 79 4 . (1) 

For ujp close to 2wo and for comparatively small F, 

5u) = \) F -u)o, \8w\<^ujq, (2) 
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even where the oscillator becomes bistable its nonlinear- 
ity remains relatively small, |7(g ,2 )| <C o> 2 . For concrete- 
ness we set 7 > 0; the results for 7 < can be obtained 
by replacing Suj — > —Su> in the final expressions. 

To describe a weakly nonlinear oscillator it is conve- 
nient to make a canonical transformation from q and p 
to the slowly varying coordinate Q and momentum P, 



U^pU 



Cpar [P cos(w F t/2) - Q sin(w F t/2)] , (3) 
-C paT u F [Psm(u F t/2) + Qcos(uj F t/2)} /2, 



where C par = (2P/37) 1 / 2 and 



[P,Q]=-iX, X = 3jH/Fuj f . 



(4) 



The dimcnsionless parameter A plays the role of H in 
quantum dynamics in the rotating frame [bj | . 

The transformed oscillator Hamiltonian has the form 
(F 2 /6 7 )5, where g = g(Q,P), 
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g(Q, P) = ^ (P 2 + Q 2 ) z + -(1 - n)P 2 -(1 + ^)Q 2 (5) 

[we use here a more conventional notation g(Q, P) in- 
stead of g(P,Q) used in Ref. The terms oc 
exp(Hnu>Ft) with n > 1 in g have been disregarded. 

The time-independent operator g is the scaled oscil- 
lator Hamiltonian in the rotating frame. Its eigenval- 
ues multiplied by F 2 /6~f give oscillator quasienergies, or 
Floquct eigenvalues. Formally, g is a Hamiltonian of an 
auxiliary stationary system with variables Q,P, and the 
eigenvalues of g give the energies of this system. The 
operator g depends on one parameter 



ji = 2lo f Suj/F . 



(6) 



For /i > —1, g(Q,P) has two minima located at P = 0, 
Q = ±(/i+ l) 1 / 2 . For < 1 the minima are separated by 
a saddle at P = Q = 0, as shown in Fig.[T] When friction 
is taken into account, the minima become stable states of 
pcriod-2 vibrations. The function g(Q, P) is symmetric 
as a consequence of the time translation symmetry: the 
change (P, Q) — > (— P, —Q) corresponds to shifting time 
in Eq. ([3]) by the modulation period r F . 

We assume the effective Planck constant A to be the 
small parameter of the theory, A <C 1. Then the low- lying 
eigenvalues of g form doublets. Splitting of the doublets 
is due to tunneling between the wells of g(Q,P). Since 
g(Q,P) = g(—Q,—P) is symmetric, the problem of level 
splitting seems to be similar to the standard problem of 
level splitting in a double- well potential [3]. As in this 
latter case, we will analyze it in the WKB approximation. 

The major distinction of the present problem comes 
from the difference between the structure of g(Q, P) and 
the Hamiltonian considered in Ref. [Til . The momentum 
P(Q; g) as given by equation g(Q, P) = g has 4 branches, 
with both real and imaginary parts in the classically for- 
bidden region of Q. This leads to new features of tunnel- 
ing and requires a modification of the method [l6| . 



We will consider splitting Sg of the two lowest eigen- 
values of g. Because of the symmetry, the corresponding 
wave functions ip±(Q) are 
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MQ)±M-Q)], 



(7) 



where ipi (Q) is the "single- well" wave function of the left 
well of g(Q, P) in Fig. [TJ It is maximal at the bottom 
of the well Qzo = — (1 + m) 1 ^ 2 an( i decays away from the 
well. To the leading order in A, the corresponding lowest 
eigenvalue of g is g min + g q , where g min = -(l + ^i) 2 /4 
is the minimum of g(Q,P) and g q = A(/i + l) 1 / 2 is the 
quantum correction. 

The wave function ipi(Q) is particularly simple for (i < 
0. In the classically forbidden region between the wells, 
IQI < \Qlo\> it has the form 



^ = C[-id P g]- 1/2 C MiSo(Q)/X], 
where Sq(Q) is given by the equation g(Q,dqSo) 



(8) 



rQ 

S (Q) = / P-(Q' 

J Qm+L q 



dQ', 



(9) 



P±{Q) = * 



sl/2 



1/2 



1 + Q 2 - ^ ± 2 (( 
L g = A/ 5 y 2 = AVV + l)-V 4 



We keep here only the contribution from the branch 
P-(Q), because P-{Q) is zero on the boundary of the 
classically forbidden range Qio + L q . For — /x » A and 
|Q/o + L q \ > \Q\ the action Sq(Q) is purely imaginary. 
The wave function ipi(Q) monotonically decays with in- 
creasing Q. 

The prefactor in the wave function @ is determined 
by the complex classical speed of the oscillator 



d P g = 2P_(Q) (Q 2 - A 7) 1/2 . 
The normalization constant C in Eq. ([8]), 
C=[( A t + l)/ 7 r] 1 /4 exp( _i / 4 ); 



(10) 



(11) 



is obtained by matching, in the range L q <C Q — Qw <C 
\Qio\- Eq. HJ to the tail of the Gaussian peak of ipi(Q), 
which is centered at Qiq. 

We are most interested in the parameter range fi ^ A 
where tunneling displays unusual behavior. For such fi 
the momentum P_ (Q) becomes complex in the range 
\Q\ < jS. This means that the decay of the wave func- 
tion is accompanied by oscillations. To correctly describe 
them we had to keep corrections oc g q in Eq. (|9|) . 

Wc first rewrite Eq. (J9j) in the form 



P-{Q) 



1 - (Q - fi) 



2 - r^Va _ 



3 9 /2 



l-lO 2 -^) 1 ^ 



(12) 
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Eq. (p~2|) applies for Q — Qiq ^> L 9 . It is seen that P-(Q) 
has two branching points inside the classically forbidden 
region. The closest to Qi is the point Qb r = —ft 1 / 2 ■ The 
WKB approximation breaks down for small Q + fti 1 / 2 . 
The wave function in this region can be shown to be pro- 
portional to Airy function Ai[-(Q + /i 1 / 2 )(2/i 1 / 2 /A 2 ) 1/3 ]. 
Therefore xjji oscillates with Q for positive Q + ft 1 / 2 . 

In contrast to the standard WKB theory of the turn- 
ing point, the prefactor in ipi contains two factors that 
experience branching at —ft 1 / 2 , see Eqs. ((HJ, (fTO]) . The 
full solution in the oscillation region can be obtained by 
going around —ft 1 / 2 in the complex plane following the 
prescription For A 2 / 3 >C Q + ft, 1 / 2 it gives 

ft ^2C\d P g\- 1 / 2 exp [-ImS (Q)/A] cos$(Q), 
$(Q) = $i(Q) + * 2 (Q). (13) 

The term Im Sq(Q) in the amplitude of the wave func- 
tion (Q2]) is determined by Eq. ©. The phase <J>(Q) has 
two terms. The term $i(Q) comes from the exponential 
factor in the WKB wave function (JSj> , 

$i(Q)=A- 1 / ReP_(Q)dQ, (14) 

where Re P-(Q) is given by Eq. (fT2|) in which we set 
(Q 2 - ft) 1 / 2 -> i(ft - g 2 ) 1/2 ; therefore Re P-(Q) > 0. It 
is simple to write $i and Im Sq(Q) in explicit form. 

The term $2(<5) in Eq. (fl"3"|) comes from the prefactor 
in 4n{Q), Eq. ©, 
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l + fi - Q< 



1/2 



7T 

4' 



(15) 
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FIG. 2: (Color online) The wave function of the ground state 
in the left well ipi (Q) in the oscillation region for A = 0.09 and 
/i = 0.5. The solid line shows explicit expressions (|13J )- (|15|l . 
the dashed line shows numerical results. Inset: ipi(Q) near its 
second zero with higher resolution. 



Decay and oscillations of the wave function described 
by Eq. (fl~3"|) are compared in Fig. [5] with the results of a 



numerical solution of the Schrodingcr equation gif) = grp. 
The left-well wave function was obtained numerically as 
a sum of the two lowest-eigenvalues solutions, cf. Eq. (J7J). 
In this calculation the basis of 120 oscillator Fock states 
was used. A good agreement between analytical and nu- 
merical results is seen already for not too small A = 0.09. 

The above solution allows us to find the tunnel split- 
ting Sg = g_ — g + of the symmetric and antisymmet- 
ric states ([7]). Following the standard approach for a 
symmetric double-well potential flH | we multiply the 
Schrodinger equations for the involved states gipi = giijji 
and gtp± = g±ip± by ip± and ip?, respectively, integrate 
over Q from — oo to and subtract the results. This gives 

Sg = -A 2 {2(1 -mM(0)$(0) 

-a 2 M(o)C(o) + $(oM'(o)]} 

or, with account taken of Eq. (fTB")) . 



Sg 
A 



16AV 2 (m + 1) 5 / 4 a/x 
(tt/x) 1 ^ 



cos [2$i(0)] 



(16) 



(M + l^ + Mln^" 1 / 2 [l + ^+l) 1 / 2 ]) 
2*i(0) = tt^A- 1 - l)/2 (m>A). 



Clearly, Sg may be positive or negative, that is, the sym- 
metric state may have a lower or higher quasiencrgy than 
the antisymmetric state. 

The dimensional splitting (F' 2 /67)|<5g| gives twice the 
matrix element of tunneling between period-2 states of 
the oscillator. This matrix element has an exponential 
factor exp(— A/X) In addition, it contains a factor 

oscillating as a function of the scaled frequency detuning 
/i/A = 6^(0;^ — 2wo)/37ft. The oscillation period is 
A(/i/A) = 4. These oscillations are shown in Fig. [3j 

The oscillations of Sg result from the wave function 
oscillations in the classically forbidden region. This can 
be seen from the analysis of ipi(Q) near the positive-Q 
boundary of the oscillation region, Q = fi 1 /' 2 . The wave 
function for Q — fi 1 / 2 3> A is a combination of the WKB 
waves with imaginary momenta P±(Q) « i[l ± (Q 2 — 
jl) 1 / 2 ]. The coefficients in this combination can be found 
in a standard way [l6| . They are determined by the phase 
$(/i 1/ ' 2 ). Only the wave with P_(Q) contributes to the 
tunneling amplitude, since P+ remains imaginary in the 
right well of g{Q,P). For ^(/i 1 / 2 ) = {An - 3)tt/4 this 
wave has zero amplitude, leading to Sg = 0. By noting 
that ^{y}/ 2 ) = 2<f>i(0) — 7r/4, we immediately obtain 
from Eq. (jTSJ) that Sg = for /i = 2nA with integer n, in 
agreement with Fig. [3] 

The occurrence of spatial oscillations of the ground 
state wave function of the scaled Hamiltonian g does not 
contradict the oscillation theorem, because g is not a sum 
of the kinetic and potential energies and is quartic in P. 
The motion described by the Hamiltonian g(Q, P) is clas- 
sically integrable. Respectively, the quantum problem is 
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reached in the range of bistability [10[ . 

In conclusion, we used the WKB approximation to 
study the wave functions of the period-2 states of a para- 
metrically modulated oscillator. We showed that these 
wave functions can display spatial oscillations in the clas- 
sically forbidden region, in the rotating frame. These 
oscillations lead to oscillations of the matrix element of 
tunneling between the period-2 states with the varying 
frequency of the modulating field. 

We are grateful to V.N. Smclyanskiy and F. Wilhclm 
for stimulating interactions. This research was supported 
in part by the NSF through grant No. PHY-0555346. 



FIG. 3: (Color online) Scaled matrix element of tunneling 
between period-2 states as a function of the scaled detuning 
of the modulation frequency from twice the oscillator eigen- 
frequency. The solid lines show explicit expression (TTO)) . the 
dashed lines show the result of numerical calculations. Inset: 
a higher-resolution plot of |<5<?|/2 vs. fx near the zero of 5g at 
fi = 6A. The data refer to A = 0.09. 



different from dynamical tunneling in classically chaotic 
systems [13, EM ; the effect we discuss has not been 
considered for such systems, to the best of our knowledge. 

The effect is also qualitatively different from photon- 
assisted/suppressed tunneling in systems with stationary 
double-well potentials: our oscillator has a single- well po- 
tential, the bistability is a consequence of resonant mod- 
ulation, and the Hamiltonian g is independent of time. 
At the same time there is a remote similarity between the 
oscillations of the tunneling matrix element for period-2 
states and for electron states in a double-well potential 
in a quantizing magnetic field [20| . However, not only is 
the physics different, but our approach is also different 
from that in Ref . I20I; in particular, it makes it possible to 
find Sg analytically. The approach can be also extended 
to a resonantly driven Duffing oscillator, where the RWA 
Hamiltonian has a structure similar to Eq. ([5]) 0, E2] • 

Tunnel splitting can be observed by preparing the sys- 
tem in one of the period-2 states and by studying inter- 
state oscillations, cf. Refs. 0, [l^. This requires that 
the tunneling rate (5ui/2fi\)\5g\ exceed ujp/iQ, where Q 
is the oscillator quality factor. The splitting sharply in- 
creases with increasing A. It will be shown separately 
that for comparatively large A (but still for \5g\ <C g q ) 
the RWA applies and relaxation remains small provided 
Sg 2 > C\/Q with C\ < 1. Our RWA numerical results 
indicate that 8g still oscillates with [i for A = 0.25 — 0.3 
and is well described by Eq. (fTu]l for fi > 2A. The local 
peak of \Sg\ for A = 0.3 and the characteristic extremum 
of cPSg/d/i 2 for A = 0.25 occur where \Sg\ « 0.01. Such 
Sg may be large enough for detecting the effect in mod- 
ulated Josephson junctions where Q = 2360 has been 
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